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Abstract 

We present an update on total and partial branching fractions and on CP asymmetries in the semi-inclusive decay 
B —» X s t + fX. Further, we summarize our results on branching fractions and CP asymmetries for semi-inclusive and 
fully-inclusive B —> X s y decays. We present the first result on the CP asymmetry difference of charged and neutral 
B —» X s y decays yielding the first constraint on the ratio of Wilson coefficients lm(Cf /Cf). 
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1. Introduction 

The decays B —» X sr jy and B —> X 5/ j{ + C are 
flavor-changing neutral current (FCNC) processes that 
are forbidden in the Standard Model (SM) at tree 
level. However, they can proceed via penguin loops 
and box diagrams. Figure [T] shows the lowest-order 
diagrams for both processes. The effective Hamil¬ 
tonian factorizes short-distance effects represented by 
perturbatively-calculable Wilson coefficients (C,) CD El 
from long-distance effects specified by four-quark oper¬ 
ators (Oj): 

//eir = ^2,V; b Vxs.dC,yr)G / . (1) 

47T 

Here, Gp is the Fermi constant, V* b and V xs d are CKM 
elements (x = u, c , t ) and fi is the renormalization 
scale. The operators have to be calculated using non- 
perturbative methods, such as the heavy quark expan¬ 
sion ID BUS ED- In B —> X s y, the dominant contribution 
arises from the magnetic dipole operator 0 1 with a top 
quark in the loop. Thus, the branching fraction depends 
on the Wilson coefficient Cf = -0.304 (NNLL) OH. 
Via operator mixing, the color-magnetic dipole opera¬ 
tor 0% contributes in higher order with Cg ff = -0.167 
(NNLL) 0 |U. In B —* X s { + {~, the weak penguin 



Figure 1: Lowest-order diagrams for B —> X s jy (top) and B 
X,/t (bottom). 


and box diagrams contribute in addition. The vector 
part is represented by operator Oq with Wilson coeffi¬ 
cient Cf = 4.211 (NNLL) 00 while the axial-vector 
part is specified by operator O m with Wilson coefficient 
Cf® = -4.103 (NNLL) 00. Again, the top quark 
in the loop yields the most dominant contribution. New 
physics adds penguin and box diagrams with new par¬ 
ticles modifying the SM values of the Wilson coeffi- 
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cients. In addition, scalar and pseudoscalar couplings 
may contribute introducing new Wilson coefficients Cs 
and Cp. Figure [2] shows examples of new physics pro¬ 
cesses involving a charged Higgs, a chargino and neu- 
tralinos 0 ua on m ej m na. These rare decays 
probe new physics at a scale of a few TeV. 

H X g>X 


b u, c, t s,d b u,c,t s,d b u,c,t s, d 

Figure 2: New physics processes with a charged Higgs bosons (left), 
a chargino plus up-type squarks (middle) and neutralinos plus down- 
type squarks (right). 


2. Study of B -► X s t + t~ 

Using a semi-inclusive approach, we have updated 
the partial and total branching fraction measurements 
of B —> X s € + t~ modes with the full BABAR data sam¬ 
ple of 471 x 10 6 BB events. We also perform the first 
measurement of direct CP asymmetry. For measuring 
partial and total branching fractions, we reconstmct 20 
exclusive final states listed in Table [T| After account¬ 
ing for K ° modes, K° s —> n Q 7t° and n° Dalitz decays, 
they represent 70% of the inclusive rate for hadronic 
masses mx s < 1.8 GeV. Using JETSET fragmenta¬ 
tion and theory predictions, we extrapolate for the miss¬ 
ing modes and those with mx, > 1.8 GeV. We im¬ 
pose requirements on the beam-energy-substituted mass 

m E s = - p*£ > 5.225 GeV and on the energy 

difference -0.1 (0.05) < A E = E* B — EqmI 2 < 0.05 for 
X s e + e~ ( X s nn ) modes where E* B and p* B are B momen¬ 
tum and B energy in the center-of-mass (CM) frame and 
Eqm is the total CM energy. We use no tagging of the B 
decay. 

To suppress e + e~ —* qq (q — u , d, s, c) events and BB 
combinatorial background, we define boosted decision 
trees (BDT) for each q 2 bin in e + e~ and pCpT separately 
(see Table 2). From these BDTs, we determine a likeli¬ 
hood ratio (Lp) to separate signal from qq and BB back¬ 
grounds. We veto J/t// and i/A2S) mass regions and use 
them as control samples. Figures[3]and [4]show the m E s 
and Lr distributions for e + e~ modes in bin <75 and for 
p' p modes in bin < 7 !, respectively. 

We measure d!B(B —> X s £ + £~)/dq 2 in six bins of 
q 2 = m 2 ( and four bins of nig defined in Tablej2] We ex¬ 
tract the signal in each bin from a two-dimensional fit to 
m E s and Lr. Figure [5] shows the differential branching 
faction as a function of q 2 (top) and m x (bottom) EE 


Table 1: Exclusive modes used in the semi-inclusive B —> X s £ + C 
analysis. 


Mode 

Mode 

B° -> K"p + p- 

B + -> K + p + p~ 

1 

+ 

O to 

T 

B + -> K + e + e~ 

B° -> K* 0 (K° s 7t°)p + p- 

B+ K ,+ (K + 7r°)ii + p- 

B° -> K*°(K + n-)p + /r 

B+ _> K* + (K° s 7T + )p + p- 

B° -> K*°(K° s 7T 0 )e + e ” 

B+ -> K* + (K + n 0 )e + e~ 

B° -> K*°(K + n-)e + e- 

B+ _> K* + (K 0 s n + )e + e- 

B° -> K° s n + 7i-)p + p- 

B + -> K° s n + 7r°p + p- 

B° -> K + n~n°p + p- 

B + -> K + n + n~p + p~ 

B° -> K° s n + 7T-)e + e~ 

B + -> K Q s n + 7T 0 e + e~ 

B° -> K + 7r-jr°e + e- 

B + — > K + n + n~e + e~ 


Table 2: Definition of the q 2 bins. 


q * 2 * bin 

q- range [GeV 2 /c 4 ] 

mu range [GeV/c 2 ] 

nix l bii 

0 

1.0 <q 2 <6.0 

1.00 < ma < 2.45 


1 

0 .1< q 2 <2.0 

0.32< met <1.41 

1 

2 

2.0< q 2 <4.3 

1.41 < m tt <2.07 

2 

3 

4.3< q 2 <8.1 

2.07 < met <2.6 

3 

4 

10.1 <q 2 <12.9 

3.18< m {{ <3.59 

4 

5 

14.2 < q 2 < ( m B - m* K ) 2 

3.77 < met < (m B - m* K ) 





Figure 3: Distributions of tues (left) and likelihood ratio (right) for 
B —» X s e + e~ in q 2 bin q$ showing data (points with error bars), the to¬ 
tal fit (thick solid blue curves), signal component (red peaking curves), 
signal cross feed (cyan/lgrey curves), BB background (magenta/dark 
grey smooth curve), e + e~ —» qq background (green/grey curves) and 
charmonium background (yellow/light grey curves). 
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Figure 4: Distributions of (left) and likelihood ratio (right) for 
B —» X s p + p~ in q 2 bin q\ showing data (points with error bars), the to¬ 
tal fit (thick solid blue curves), signal component (red peaking curves), 
signal cross feed (cyan/lgrey curves), BB background (magenta/dark 
grey smooth curve), e + e~ —» qq background (green/grey curves) and 
charmonium background (yellow/light grey curves) 


Table 3: The B —> X s l + £ branching fraction measurements in the 
low and high q 2 regions m in comparison to the SM prediction. 


Mode 

BABAR [10 b ] 

SM [10^ b ] 

^ 2 [GeV 2 /c 4 ] 

1-6 

1-6 

B —> X s n + iT 

0 66 +0 - 82+0 - 30 + 0 07 

-0.76-0.24 - u,u/ 

1.59 + 0.11 

B -+ X s e + e- 

2 93+0.47+0.21 + oi8 

-0.45-0.16 - u ' io 

1.64 + 0.11 

b -> x s e + r 

1 60 +a41+017 + 0 07 

1,u -0.39-0.13 


r[GeV 2 /c 4 ] 

> 14.2 

> 14.2 

B -+ X s v + R- 

0 60 +0 - 31+0 05 + 0 00 
-0.29-0.04 - u,uu 

0 25 +0 07 
u ‘ z ' -0.06 

B -> X s e + e~ 

0 56 +0 - 19+0 03 + 0 00 

U-JU -0.18-0.03 - u - uu 


B -> x s { + r 

0 57 +016+0 - 03 + 0 00 

7 -0.15-0.02 - u,uu 



Table [3] summarizes the differential branching fractions 
in the low and high q 2 regions in comparison to the SM 
predictions lfl7l fl8l fl9l fl8l l20l I2T11^21 l23l l24l 1251 l26l 
|27|. In both regions of q 2 , the differential branching 
fraction is in good agreement with the SM prediction. 
These results supersede the previous BABAR measure¬ 
ments ll28l and are in good agreement with the Belle 

results (291- 




Figure 5: Differential branching fraction of B —> X s e + e~ (blue 
points), B —> X s q' q (black squares), and B —> X S C*C~ (red trian¬ 
gles) versus q 2 (top) and versus mx s (bottom) in comparison to the 
SM prediction (histogram). The grey-shaded bands show the J/ip and 
i//(2S) vetoed regions. 


The direct CP asymmetry is defined by: 

s (b -> x s ee~) - s(b -> x s ez~) 

fficp — —=---• ( 2 ) 

£{B -4 X s { + i-) + £{B -> X s e + e~) 

We use 14 self-tagging modes consisting of all B + 
modes and the B° modes with decays to a K' listed in 
Table [T] to measure C p(B —* X s £ + t~) in five q 1 bins. 
Note that we have combined bins <74 and <75 due to low 
statistics. Figure [ 6 ] shows the CP asymmetry as a func¬ 
tion of q 1 . The SM prediction of the CP asymmetry in 
the entire q 2 region is close to zero (30) [31] [32] 0 ■ In 
new physics models, however, cp may be significantly 
enhanced ItTTl [33ll . In the full range of q 2 we measure 
$[ C p = 0.04 + 0.11 +0.01 |T6l , which is in good agree¬ 
ment with the SM prediction. The CP asymmetries in 
the five q 1 bins are also consistent with zero. 


0.5- 




- 1 - 

0 


J/\l/ 


20 

q 2 (GeV/c 2 )' 


25 

2x2 


Figure 6: The CP asymmetry as a function of q 2 . The grey-shaded 
bands show the J/i{/ and il/(2S) vetoed regions. 
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3. Study of B —> X s y 


In the SM, the B — * X s y branching fraction is calcu¬ 
lated in next-to-next leading order (4 loops) yielding 

S(B -> X s y) = (3.15 ± 0.23) x 10^ 4 (3) 

for photon energies E y > 1.6 GeV 1: 341135 1. 

To extract the B —» X s y signal experimentally from 
e + e —> BB and e + e~ —» qq backgrounds, we use two 
very different strategies. The first strategy consists of a 
semi-inclusive approach in which we sum over 38 ex¬ 
clusive B —> X s y final states with lAT*^ lAfJ?) or 3 AT*, 
< 47t(< 27t°), and < 1?/. We use no tagging of the other 
B meson. We need to model the missing modes. Due 
to large backgrounds, we select events with a minimum 
photon energy of E y > 1.9 GeV and then extrapolate the 
branching fraction to photon energies E y > 1.6 GeV. 
With this approach, we measure the branching fraction, 
CP asymmetry and the difference in CP asymmetries be¬ 
tween charged and neutral B decays using 471 x 10 6 BB 
events (361 . 

The second strategy is a fully inclusive approach. To 
suppress backgrounds from BB and qq decays, we im¬ 
pose stringent constraints on isolated photons to remove 
clusters that may have originated from 7r° and i] decays. 
We use a semileptonic tag of the other B meson and re¬ 
quire a minimum photon energy of E y > 1.8 GeV but 
impose no requirements on the hadronic mass system. 
Using 383 x 10 6 BB events, we measure the B —» X s y 
branching fraction measurement and the CP asymmetry 

for B -» X s+d y EZJED- 

Table [4] summarizes our B —» X s y branching frac¬ 
tion measurements of the semi-inclusive and fully inclu¬ 
sive methods (361 E3 [38] . Figure [7] shows the BA BAR 
results extrapolated to a minimum photon energy of 
1.6 GeV in comparison to results from Belle Ii40l 14T1 
1421 . CLEO 03] and the SM prediction [gj|35]. Our 
results are in good agreement with those of the other 
experiments as well as the SM prediction. 

For the semi-inclusive method, the direct CP asym¬ 
metry is defined by: 


/A cp(X s y) = 


B(B -> X s y) - B(B -> X s y) 
B(B X s y) + B(B ^ X s y)' 


(4) 


The SM prediction yields -0.6% < /Acp(B —> X s y) < 
2.8% (321461. Using 16 self-tagging exclusive modes 
and 471 x 10 6 BB events, we measure /A CP (B —> X s y) = 
(1.7 ± 1.9** ± 1.0 JVJ )% B2- This supersedes the old 
BAJ3AR measurement l48l .We further measures the CP 
asymmetry difference between charged and neutral B 
decays: 

A # C p = & CP (B + -+ X+y) - Wcp(B° -» X° s y), (5) 


Table 4: Our measurements of !B{B —> X s y ) from the semi- 

inclusive EU and fully-inclusive ED analyses and their extrapola¬ 
tions to Ey > 1.6 GeV. The first uncertainty is statistical, the second 
is systematic and the third is from model dependence and extrapola¬ 
tion to 1.6 GeV. 


method 

E y > 

B(B -> X s y) [10- 4 ] 

semi¬ 

exclusive 

1.9 GeV 
1.6 GeV 

3.29 + 0.19 + 0.48 

3.52 + 0.20 + 0.51 +0.04 

inclusive 

1.8 GeV 
1.6 GeV 

3.21+0.15 + 0.29 + 0.08 
3.31+0.16 + 0.30 + 0.10 



Figure 7: Summary of !B(B —> X s y) measurements from BABAR El 
E3EI1E11. Belle mmm, CLEO El and the HFAG average El 
in comparison to the SM prediction EHED after extrapolation to 
E*> 1.6 GeV. 


which depends on the Wilson coefficients C® ff and CJj IT : 


A78 

m h 


c: 


' e ff r 

A/Acp = 4 ir-as—Im-^ - 0.12--— Im — a : (6) 

LF Cf 100 MeV cf 


where the scale parameter A 7 g is constrained by 
17 MeV < A 78 < 190 MeV. In the SM, Cf and Cf are 
real so that A/Acp vanishes. However in new physics 
models, these Wilson coefficients may have imaginary 
parts yielding a non-vanishing A/A cp- 

From a simultaneous fit to charged and neutral B de¬ 
cays, we measure A/Acp(B — » X s y) = (5.0 ± 3.9 stat ± 
1.5 jys )% from which we set an upper and lower limit at 
90% CL on Im(Cf/Cf) 03: 

£ e ff 

- 1.64 < Im-^= < 6.52 at 90% CL. (7) 

/-<eff 

This is the first A/Acp measurements and the first con¬ 
straint on 7m(Cg ff /C® ff ). Figure[8](top) shows the Ay 1 of 
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the fit as a function of ImiCf /Of). The shape of Ay 2 
as a function of Jm(Co ff /CS ff ) is not parabolic indicat¬ 
ing that the likelihood has a non-Gaussian shape. The 
reason is that Ax 2 is determined from all possible val¬ 
ues of A 78 . In the region ~ 0.2 < <~ 

2.6 a change in J/n(Cg ff /CS ff ) Ax 1 can be compen¬ 
sated by a change in A 78 leaving Ax 1 unchanged. For 
positive values larger (smaller) than 2.6 (0.2), Ax 1 in¬ 
creases slowly (rapidly), since A 7X remains nearly con¬ 
stant at the minimum value (increases rapidly). For neg¬ 
ative Jw(Cg ff /CS ff ) values, A 78 starts to decrease again, 
which leads to a change in the Ax 2 shape. Figure[8](bot- 
tom) shows A 78 as a function of 7w(Cg ff /C® ff ). 



>< 



Im 7T 


Figure 8: The A\ 2 function versus Jm(Cg ff /C^ ff ) (top) and the de¬ 
pendence of A 78 on Im(C < f /C c 1 a ) (bottom). The blue dark-shaded 
(orange light-shaded) regions show the 68% (90%) CL intervals. 


In the fully-inclusive analysis, the B —* X d decay can¬ 
not be separated from the B —> X s decay and we mea¬ 
sure: 


■ftcp(x s+d y) - 


S(B 

Wfi 


X s+d y) - S(B -> X s+d y) 
X s+d y) + S(B ^ X s + d y) 


In the SM, ABcpiB —> X s+d y) is zero 11491 . From the 
charge of the B and B, we determine the CP asymme¬ 
try. Using 383 x 10 6 BB events, we measure £Bcp(B —> 
X-s+dy) = (5-7 ± 6.0 ± 1.8)%, which is consistent with 
the SM prediction J49[. Figure [9] shows a summary of 
all CP asymmetry measurements in comparison to the 
SM predictions. 


SM A(B-»X s+d y) prl 106,141801 (2011) 



Figure 9: Summary of cp measurements for B —> X s y from semi- 
inclusive analyses (BABAR 1471 , Belle M501I ) and for B —> X s+ dy 
from fully inclusive analyses (BABAR | 37] |38] [39], CLEO PTIl ), 
Belle (52) and the HFAG average ED in comparison to the SM pre¬ 
diction for B —> X s y (mono. 


4. Conclusion 


We performed the first cp measurement in five q 1 
bins in semi-inclusive B —» X s £ + (~ decays and up¬ 
dated the differential branching fraction. The B —> 
X s £ + C~ partial branching fractions and CP asymmetries 
are in good agreement with the SM predictions. Our 
lB(p measurement in the semi-inclusive B —» X s y de¬ 
cay is the most precise CP asymmetry measurement. 
The A£Bcp(B —> X s y) result yields first constraint on 
7;w(Cg/C®). The B —> X s y branching fractions and 
CP asymmetries are both in good agreement with the 
SM predictions. New progress on these inclusive de¬ 
cays will come from Belle II. For the B —» X s y and 
B —» X s l*C semi-inclusive decays, we expect preci¬ 
sion measurements. For the inclusive B —> X s y and 
B —» X s £ + { decays, we expect new possibilities by tag¬ 
ging the other B meson via full B reconstruction. 
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